Method for segmenting images with intensity-based label propagation and probabilistic of level sets

ABSTRACT

A method exploits user labels in image segmentation. The user labels are propagated with respect to image intensity information. Propagated user labels are included in a cost function of level set evolution. The level set represents a probability of the object segment.

FIELD OF THE INVENTION

This invention relates generally to image processing, and more particularly to image segmentation.

BACKGROUND OF THE INVENTION

Conventional image processing can use a level set and user information to segment an image, e.g., see U.S. Pat. No. 7,277,582, “User interactive level set methods for image segmentation.” The level set is evolved to minimize a cost function, wherein variances of image intensities object pixels and background pixels are used in cost function. The user can label pixels as object or background pixels and penalize deviations of the sign of the evolving level set from the sign of the user-labeled pixels, and nearby pixels with an isotropic weighting, i.e., the weighting is uniform in all directions.

However, conventionally, the user labels do not incorporate image intensity information. The variance terms in the level set evolution dominate and only pixels near user selected background locations remain as background.

FIG. 1 shows an image that has a failure of the segmentation with conventional user interaction. The correct segmentation shown in FIG. 2 is difficult to achieve.

It is an object of the invention to produce a segmentation that approximated the correct segmentation.

SUMMARY OF THE INVENTION

The embodiments of the invention provide a method for segmenting an input image using level sets and user labels. The method propagates the labels based on image intensity information, and combines the information with a probabilistic level set representation. The method improves increases accuracy, while reducing processing time.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is an image of a failed segmentation with conventional methods;

FIG. 2 is the correct segmentation for the image in FIG. 1;

FIG. 3 is a block diagram of a method for segmenting an image according to embodiments of the invention;

FIG. 4 is an image with labeled pixels according to embodiments of the invention; and

FIG. 5 is a schematic of an image grid graph according to embodiments of the invention; and

FIG. 6 is an image of a segment according to embodiments of the invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

FIG. 3 shows a method 300 for segmenting an input image 301 according to embodiments of our invention to produces a segment 309 in an output image 309, see FIG. 6. The steps of the method can be performed in a processor including memory and input/output interfaces connected by a bus as known in the art.

A user supplies a set of labels 310. The set of labels identifies object pixels 311 and background pixels 312, see FIG. 4. The labels are propagated 310 to produce an initial level set 303. Then, the level set is evolved, using a cost function J described in detail below, to produce a segment 309 in and output image.

The level set represents the segmentation of the input image in a spatial domain with functions Ω₁, Ω₂ and C, where Ω₁={(x,y): ƒ(x,y)>0}, Ω₂={(x,y): ƒ(x,y)<0}, C={(x,y): ƒ(x,y)=0}, where the function Ω₁ identifies a background region of the object pixels, the function Ω₂ identifies a background region of the background pixels, and the function C identifies boundary pixels between the object and the background regions, and (x, y) are coordinates of locations of the pixels in the input image.

A region-based level set cost function J_(o) of the cost function J is

$\begin{matrix} {{J_{0} = {\sum\limits_{i = 1}^{2}{\int_{\Omega}{\left( {{u\left( {x,y} \right)} - {\overset{\_}{u}}_{i}} \right)^{2}{\chi_{i}\left( {f\left( {x,y} \right)} \right)}{\mathbb{d}x}{\mathbb{d}y}}}}},} & (4) \end{matrix}$

where u is an intensity of the pixel at image location (x, y), and

$\begin{matrix} {{{\overset{\_}{u}}_{i} = \frac{\int_{\Omega}{{u\left( {x,y} \right)}\chi_{i}\;\left( {f\left( {x,y} \right)} \right){\mathbb{d}x}{\mathbb{d}y}}}{\int_{\Omega}{{\chi_{i}\left( {f\left( {x,y} \right)} \right)}{\mathbb{d}x}{\mathbb{d}y}}}},} & (5) \end{matrix}$ are means of the pixel intensities of each region, and

$\begin{matrix} {{{\chi_{1}(z)} = {H(z)}},{{\chi_{2}(z)} = {1 - {H(z)}}},{where}} & (6) \\ {{H(z)} = \left\{ \begin{matrix} 0 & {{{if}\mspace{14mu} z} < 0} \\ 1 & {{{{if}\mspace{14mu} z} \geq 0},} \end{matrix} \right.} & (7) \end{matrix}$ is a Heaviside step function. The Heaviside step function H, also called a unit step function, is a discontinuous function whose value is zero for negative argument, and one for positive arguments.

Our cost function J also includes a user label cost function J_(U) and a maximum belief cost function J_(B) as described below.

The set of user labels is encoded using a function

$\begin{matrix} {{L\left( {x,y} \right)} = \left\{ \begin{matrix} 1 & \begin{matrix} {{if}\mspace{14mu}{the}\mspace{14mu}{pixel}\mspace{14mu}{at}\mspace{14mu}{the}\mspace{14mu}{location}\mspace{14mu}\left( {x,y} \right)} \\ {{{is}\mspace{14mu}{one}\mspace{14mu}{of}\mspace{14mu}{the}\mspace{14mu}{object}\mspace{14mu}{pixels}},} \end{matrix} \\ 0 & \begin{matrix} {{if}\mspace{14mu}{the}\mspace{14mu}{pixel}\mspace{14mu}{at}\mspace{14mu}{the}\mspace{14mu}{location}\mspace{14mu}\left( {x,y} \right)} \\ {{{is}\mspace{14mu}{one}\mspace{14mu}{of}\mspace{14mu}{the}\mspace{14mu}{background}\mspace{14mu}{pixels}},{and}} \end{matrix} \\ 0.5 & {{otherwise}.} \end{matrix} \right.} & (8) \end{matrix}$

As shown in FIG. 5, the label propagation 310 determines a minimum distance D along a path 501 in an image grid graph from each unlabeled pixel 510 to one of the labeled pixel 511. In the image grid graph, each pixel is a node and edges connect each node with eight nearest neighboring pixels as shown only for pixel 510.

The minimum distance is

$\begin{matrix} {{{D_{k}\left( {x,y} \right)} = {\min\limits_{{({x_{1},y_{1}})},{({x_{2},y_{2}})},\mspace{14mu}\ldots\mspace{14mu},{({x_{m},y_{m}})}}{\sum\limits_{i = 1}^{m}{d\left( {\left( {x_{i - 1},y_{i - 1}} \right),\left( {x_{i},y_{i}} \right)} \right)}}}},} & (9) \end{matrix}$ where L(x_(m), y_(m))=k, for k=0, 1, and where m is an index of the pixel (x_(m), y_(m)).

That is, D₀(x, y) is the minimum distance from the unlabeled pixel at location (x, y) to one of the labeled pixels that is similar in intensity. The corresponding path is (x, y),(x ₁ , y ₁), (x ₂ , y ₂), . . . , (x _(m) , y _(m)), where L(x_(m), y_(m))=0. The index m changes for every pixel location (x, y).

In other words, unlabeled pixels at different locations have different paths to one of the labeled pixel, and hence different path lengths m. d((X _(i−1) , Y _(i−1)),(X _(i) , Y _(i)))=∥U(X _(i−1) , Y _(i−1))−U(X _(i) , Y _(i))∥  (10)

The minimum distances initialize a probabilistic level set as

$\begin{matrix} {{{f^{0}\left( {x,y} \right)} = {1 - \frac{D_{1}\left( {x,y} \right)}{{D_{1}\left( {x,y} \right)} + {D_{0}\left( {x,y} \right)}}}},{and}} & (11) \\ {J_{U} = {- {\int_{\Omega}{\left( {{f\left( {x,y} \right)} - {f^{0}\left( {x,y} \right)}} \right)^{2}{\mathbb{d}x}{{\mathbb{d}y}.}}}}} & (12) \end{matrix}$

The probabilistic level sets is used in the maximum belief cost function:

$\begin{matrix} {{J_{B} = \frac{\int_{\Omega}{{f\left( {x,y} \right)}{H\left( {f\left( {x,y} \right)} \right)}{\mathbb{d}x}{\mathbb{d}y}}}{N}},} & (13) \end{matrix}$

where N=∫ _(Ω) H(ƒ(x, y))dxdy is a number of the object pixels.

Thus, the cost function used during the evolution is the sum of the region-based cost function J_(o), the user cost function J_(U), and the maximum belief cost function J_(B): J=J _(O)+λ_(U) J _(U)+λ_(B) J _(B)  (14) where λ_(U)_and λ_(B) are scalar parameters to control a relative importance of the user cost function J_(U) and the maximum belief cost function J_(B), respectively.

The level set is evolved using Euler-Lagrange equations of the cost function until a local minimum is reached. The EulerLagrange equations are differentiable and stationary at the local minimum.

FIG. 6 shows the segment 309 produce by the method described above, which is visible more accurate than the segmentation shown in FIG. 1 with conventional methods.

Although the invention has been described by way of examples of preferred embodiments, it is to be understood that various other adaptations and modifications may be made within the spirit and scope of the invention. Therefore, it is the object of the appended claims to cover all such variations and modifications as come within the true spirit and scope of the invention. 

1. A method for segmenting an input image of an object, comprising the steps of: specifying, by a user, a set of labels for pixels in the input image to produce labeled pixels and remaining unlabeled pixel, wherein the set of labels identifies object pixels and background pixels; propagating the labels in the input image based on differences of intensities of neighboring pixels to produce an initial level set, wherein the initial level set is represented probabilistically; and evolving the initial level set using a cost function to produce a segment in an output image, and wherein the steps are performed in a processor.
 2. The method of claim 1, wherein the level represents a segmentation of the input image in a spatial domain with functions Ω₁, Ω₂ and C, where Ω₁={(x,y): ƒ(x,y)>0}, Ω₂={(x,y): ƒ(x,y)<0}, C={(x,y): ƒ(x,y)=0}, and where the function Ω₁ identifies a object region of the object pixels, the function Ω2 identifies a background region of the background pixel, and the function C identifies boundary pixels between the object regions and the background region, and (x, y) are coordinates of locations of the pixels in the input image.
 3. The method of claim 1, wherein the cost function includes a region-based level set cost function J_(o), a user label cost function J_(U), and a maximum belief cost function J_(B).
 4. The method of claim 3, wherein the region-based cost function is ${J_{0} = {\sum\limits_{i = 1}^{2}{\int_{\Omega}{\left( {{u\left( {x,y} \right)} - {\overset{\_}{u}}_{i}} \right)^{2}{\chi_{i}\left( {f\left( {x,y} \right)} \right)}{\mathbb{d}x}{\mathbb{d}y}}}}},$ where u is the intensity of the pixel at the location (x, y), and $\begin{matrix} {{{\overset{\_}{u}}_{i} = \frac{\int_{\Omega}{{u\left( {x,y} \right)}{\chi_{i}\left( {f\left( {x,y} \right)} \right)}{\mathbb{d}x}{\mathbb{d}y}}}{\int_{\Omega}{{\chi_{i}\left( {f\left( {x,y} \right)} \right)}{\mathbb{d}x}{\mathbb{d}y}}}},} & (5) \end{matrix}$ are means of the pixel intensities in a particular region, and χ₁(z) = H(z), χ₂(z) = 1 − H(z), where ${H(z)} = \left\{ \begin{matrix} 0 & {{{if}\mspace{14mu} z} < 0} \\ 1 & {{{if}\mspace{14mu} z} \geq 0} \end{matrix} \right.$ is a Heaviside step function.
 5. The method of claim 4, where the set of user labels is encoded using a function ${L\left( {x,y} \right)} = \left\{ \begin{matrix} 1 & \begin{matrix} {{if}\mspace{14mu}{the}\mspace{14mu}{pixel}\mspace{14mu}{at}\mspace{14mu}{the}\mspace{14mu}{location}\mspace{14mu}\left( {x,y} \right)} \\ {{{is}\mspace{14mu}{one}\mspace{14mu}{of}\mspace{14mu}{the}\mspace{14mu}{object}\mspace{14mu}{pixels}},} \end{matrix} \\ 0 & \begin{matrix} {{if}\mspace{14mu}{the}\mspace{14mu}{pixel}\mspace{14mu}{at}\mspace{14mu}{the}\mspace{14mu}{location}\mspace{14mu}\left( {x,y} \right)} \\ {{{is}\mspace{14mu}{one}\mspace{14mu}{of}\mspace{14mu}{the}\mspace{14mu}{background}\mspace{14mu}{pixels}},{and}} \end{matrix} \\ 0.5 & {{otherwise}.} \end{matrix} \right.$
 6. The method of claim 5, wherein the propagating determines a minimum distance D along a path in an image grid graph from each unlabeled pixel to one of the labeled pixel, wherein the minimum distance is $\begin{matrix} {{{D_{k}\left( {x,y} \right)} = {\min\limits_{{({x_{1},y_{1}})},{({x_{2},y_{2}})},\mspace{14mu}\ldots\mspace{14mu},{({x_{m},y_{m}})}}{\sum\limits_{i = 1}^{m}{d\left( {\left( {x_{i - 1},y_{i - 1}} \right),\left( {x_{i},y_{i}} \right)} \right)}}}},} & (9) \end{matrix}$ where L(X_(m), Y_(m))=k, for k=0, 1, and m is an index of the pixel (X_(m), Y_(m)).
 7. The method of claim 6, wherein the minimum distances initialize a probabilistic level set as ${{f^{0}\left( {x,y} \right)} = {1 - \frac{D_{1}\left( {x,y} \right)}{{D_{1}\left( {x,y} \right)} + {D_{0}\left( {x,y} \right)}}}},{and}$ J_(U) = −∫_(Ω)(f(x, y) − f⁰(x, y))²𝕕x𝕕y.
 8. The method of claim 7, wherein the probabilistic level sets is used in the maximum belief cost function: ${J_{B} = \frac{\int_{\Omega}{{f\left( {x,y} \right)}{H\left( {f\left( {x,y} \right)} \right)}{\mathbb{d}x}{\mathbb{d}y}}}{N}},$ where N=∫ _(Ω) H(ƒ(x,y))dxdy is a number of the object pixels.
 9. The method of claim 1, wherein the cost function is J=J _(O)+λ_(U) J _(U)+λ_(B) J _(B)  (14) where λ_(U) and λ_(B) are scalar parameters to control a relative importance of the user cost function J_(U) and the maximum belief cost function J_(B), respectively.
 10. The method of claim 1, wherein the evolving uses Euler-Lagrange equations of the cost function until a local minimum is reached. 